Magnetic field induced band insulator to Mott insulator transformations in 
4-component alkali fermions at half-filling 
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Under the influence of an external magnetic field and spin-changing collisions, the band insulator 
(BI) state of one-dimensional (ID) s-wave repulsively interacting 4-component fermions at half-filling 
transforms into Mott insulator (MI) states with spontaneously broken translational symmetry: a 
dimerized state for shallow lattices and a Neel state for deep lattices via an intermediate topological 
state. Since a BI has vanishing entropy per particle, these MI phases could be particularly inviting 
for experimental realization under the similar conditions as those for 40 K atoms [J, provided the 
magnetic field is changed adiabatically. 



Ultra-cold spinor Fermi gases in optical lattices are 
a potentially extraordinary tool for the modelling of 
strongly-correlated condensed matter systems under 
well-controlled conditions [U, Q . The metal-to- MI tran- 
sition was experimentally observed with ultra-cold 2- 
component Fermi gases [l, H[. The MI phase should 
acquire a magnetic Neel (antifcrromagnctic) ordering at 
sufficiently low temperatures, due to the super-exchange 
interaction, which at the moment has not yet been re- 
solved experimentally [f|, the main obstacle being the 
absence of efficient cooling methods in strongly corre- 
lated regimes in the presence of an optical lattice Q- 
With increasing number of components new interesting 
physics emerges as multi-component lattice Fermi gases 
offer possibility of hosting exotic ground state, like 2D 
spin liquids [|| or ID topological states fol, Hoj|. 

Dynamical properties of multi-component ultra-cold 
fermions have attracted recent experimental and theo- 
retical attention as well. Spinor dynamics were studied 
in the case of s-wave interacting 4-component fermions 
both experimentally llf and theoretically within mean- 
field approximation [ill ]. In recent experiment [l| ini- 
tially a state of 2 atoms per site was stabilized by the 
quadratic Zeeman effect; the starting state being a BI of 
the 2-component Fermi Hubbard model. Although the 
40 K atoms were loaded in a 3D lattice, hopping was al- 
lowed along one direction, and since s-wave interactions 
arc only felt on-site the system studied experimentally 
was collection of ID decoupled lattices. After quench- 
ing the magnetic field, multiflavour spin dynamics and 
strong damping due to many-body effects were observed. 

In this letter we analyse the ground state phases of 4- 
component alkali fermions at half-filling using non- per- 
turbative analytical and numerical tools, particularly tai- 
lored for studying ID systems. In the parameter space of 
magnetic field and lattice depth we identify various MI 
phases: those, that break translational symmetry and 
are characterised by local order and a topological Hal- 
dane phase which does not break any microscopic symme- 
try and is characterised by non-local string order. Even 
though these different MI phases arc manifested below 
the ultra-low temperatures the fact that one can start 
from the BI state which has vanishing entropy per par- 
ticle, and by adiabatically changing the magnetic field 
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FIG. 1: (Color online) Numerical ground state phase dia- 
gram for a 4-component mixture of 40 K atoms at half-filling 
with average interaction energy U ~ 4 ^ (a^ + + af + 
+ ag")/5 (a^,...a^ being the s-wave scattering lengths 
of 40 K) and quadratic Zeeman coupling q (both in units of 
hopping t set by lattice depth). Dotted line indicates sharp 
cross-over and continuous lines are quantum phase transition 
lines, crossing at multicritical point M. For nature of different 
phases and phase transitions consult the text. 



enter into these non-trivial MI phases, makes the system 
of half-filled 4-componcnt alkali fermions a particularly 
attractive candidate to resolve ground state spin order 
in experiments on ultra-cold lattice gases. Note that the 
spin-changing collision processes, crucial for BI to MI 
transformations, are maximally pronounced at half-filling 
corresponding to occupancy of 2 particles per site. 

Our main result, the ground state phase diagram of 
4 components of 40 K atoms, is presented in FigJT] The 
most relevant finding is that for deep lattices a topolog- 
ical Haldane state developes next to the BI state. This 
can be important for experiments because one can ar- 
rive at a topological state starting from the BI of the 
two component mixture and activate spin-changing colli- 
sions by adiabatically changing the magnetic field. This 
phase diagram is not particular of 40 K atoms, it captures 
generic phases of repulsive 4-component alkali fermions 
at half-filling, provided these components do not form 
complete hyperfine multiplet. For the case of hyperfine 
spin / = 3/2, due to large internal symmetry which holds 
even for nonzero quadratic Zeeman coupling [l2| , the sit- 
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uation simplifies and only dimer and singlct/BI phases 
are realized as depicted in FigOJb). 

f=3/2 multiylet- We start from the case of a hypcrfmc 
spin / = 3/2 The lattice Hamiltonian that describes 
s-wave interacting fermions in the presence of quadratic 
Zeeman coupling q is a sum of parts quadratic and quartic 
in fermions, 

L L 

h = H ,j + H in t,j, (i) 

3/2 
a=-3/2 

HintJ = 9F PF, mF (J) P F,m F (J)- 

F,\m F \<F 

The operator ip a annihilates an atom in the hyper- 
fine spin state \ f = 3/2, rrif = a) and n a ,j — '4 ! aj'4'a,j- 
Interaction coefficients qp are proportional to scattering 
strengths ~ 4irh 2 a,F /m [141 ]. where m is the atomic mass, 
clf are the s-wave scattering lengths and Pp m (j) = 

^2 a p(F,fnF\a0)ipl l jip'g ■ create on-site pairs with total 
spin F and projection mp. 

Due to Pauli principle only channels with F = 0, 2 
are allowed in the low energy scattering. Interactions are 
assumed to be repulsive and we concentrate on the case 
of two fermions per lattice site on average and zero net 
magnetization N a = N— a , where N a — Ylj( n ct,j}- 

Usually, in alkali metals, differences between the s- 
wave scattering lengths are a few percent of the av- 
erage scattering length, hence deviations from SU(4) 
symmetry (ao = a-2 case) are small. At half- 
filling a unique gapped phase, with spontaneously dou- 
bled lattice constant and dimer (spin-Peierls) order, 
is realized with an order parameter defined as D = 
IE Qj (-l) j (V^VW+i)!/ (L - 1) > 0, in the vicinity of 
-St/ (4) point in the absence of magnetic field [lf|[l6|. We 
show that for ao ^ a,i the magnetic field above a critical 
value restores translational symmetry via quantum phase 
transition belonging to the second order Ising universal- 
ity class and the system enters a site-singlet state that is 
adiabatically connected to the BI of the two-component 
Fermi Hubbard model. For the case ao = a2, relevant for 
alkaline-earth atoms [13], the dimer state is separated 
from BI with gapless Luttinger liquid state. 

For the 2-component repulsive Hubbard model at half 
filling the charge sector is gapped and low energy effec- 
tive theory is described by spin-1/2 Heiscnbcrg chain. In 
the case of 4-component fermions, at half filling, sim- 
ilarly the charge sector is gapped and instead of one 
spin sector there are 3 flavour sectors described by the 
dual bosonic fields [<f> v , 6 V ], [fai, 9 t i], and [<p t 2, 0a], 
with [0(x),d y (jj\_= i5{x — y), all of which are generi- 
cally gapped The quadratic Zeeman effect cou- 

ples only to the chiral sector cj> v . Chirality is defined 

as T = ^Ej([ n |,j + n -lj ~ n U ~ n -§,j}) ~ i d( t>v), 



and it measures imbalance between ±i and ±| compo- 
nents. At half-filling with a relevant 4kp Umklapp, the 
important low energy sector is governed by the following 
Hamiltonian density, 

n v = y |^(W 2 + -^(9 X ^) 2 | 

+ g cos V^TT(f>v + g-y cos VAtt9 v — d x (j> v , (2) 



where, v v = ^vp + (g - 3g 2 )vF K v = 

^/2ttvf / {2ttvf + go — 332), vf being the Fermi velocity 
of free fermions. Coupling constants g ~ go + <?2 and 
<?7 ~ 52 — So 7 include proportionality coefficients involv- 
ing averages of massive charge and <fit 1:2 fields [i"8j . 

The effective model in Eq.([2]) is diagonalized by 
refermionization procedure [l9j and Bogoliubov transfor- 
mation. The low energy excitation spectrum allows us to 
locate quantum phase transition and determine its na- 
ture. In general for g 7 7^ the excitation spectrum is 
gapped, except at a critical value of the q = q c . The low 
energy spectrum is given by two fermionic branches, 

w±(A) = \Jg 2 +g? l + V + fc 2 ± 2^ 2 g 2 + 4o 2 ( 5 2 + &). 

Only one branch becomes gapless for q c = g 2 — g 2 /2, 

cj_(fc) = V-\k\ + 0(k 2 ). The velocity of the linear mode, 
at the criticality is u_ = g-y/g and the transition belongs 
to the 2D Ising universality class [20j| . For g 7 = (with 
additional t/(l) symmetry corresponding to conservation 
of chirality), at q = q c = g/2, the linear velocity vanishes 
and the Ising transition transforms into commensurate- 
incommensurate one [l9j with characteristic quadratic 
(non-relativistic) low energy dispersion ui_(k) ~ k 2 /2g. 
Dimer phase that is realized for q < q c transforms into 
singlet phases for q > q c . Two singlet phases, corre- 
sponding to g~{ > and g-y < cases respectively, are 
separated by a Gaussian phase transition between each 
other at smaller values of q > q c , however both of them 
evolve adiabatically into BI state with increasing q. The 
singlet phases can be approximated as \S) ~ Yij 

where \S), = |±), <g> | - ±),- - C±l§>j ® I - |>j, with 
Q± — > ±1 for q —> and |C±| — > after the crossover into 
BI. The dependence of chirality and chirality susceptibil- 
ity obtained from the effective model Eq. @ at the dimer 
to singlet phase transitions for different values of g 7 is 
depicted in Fig{2j 

To check the phase diagram at strong couplings (rel- 
evant for experiments) we use density matrix renormal- 
isation group (DMRG) calculations for open boundary 
conditions [21|. We keep in average up to 800 states and 
L = 60 sites, comparable to the number of occupied lat- 
tice sites per chain We monitor directly dimer order 
parameter D. The scaling of dimer order and a collapse 
of the data for different system sizes on a single curve de- 
picted in Fig. (J3J) (a) confirms Ising character of the dimer 
to singlet phase transitions D ~ \q — g c | 1//8 . The numer- 
ical ground state phase diagram for / = 3/2 fermions at 
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FIG. 2: (Color online) Chirality versus quadratic Zeeman cou- 
pling (dashed lines) for different values of g~, (g-, = for the 
right most curve and g-f = g for the left most curve). Inter- 
sections of the continuous line with the dashed lines give q c 
and critical value of chirality. Inset shows chirality suscepti- 
bility and for right most curve the divergence corresponds to 
commensurate-incommensurate phase transition. 




0.2 q 0.3 



y 0.5 



FIG. 3: (Color online) (a) Scaling of the dimer order parame- 
ter for different system sizes. The inset shows the collapse of 
all finite-system-size data to a single curve, confirming Ising 
criticality. (b) Ground state phase diagram for / = 3/2 alkali 
fcrmions at half filling, as function of q (in units of t) and 
7 = {92 ~ ffo)/2t for the case of (g2 + 3o)/2 = Wt. 

half-filling is depicted in FigJ^b). Fidelity susceptibil- 
ity (not shown) has a pronounced peak, both at dimer 
to singlet phase transition as well as at the singlet to 
BI crossover, however for the latter case the height does 
not scale with system size, indicating absence of phase 
transition (but rather the presence of a sharp crossover) 
between singlet and BI states. In the BI state chirality 
(quasi) saturates t ~ 1. 

The Ising phase transition between dimer and singlet 
phases, induced by an external magnetic field, general- 
izes the one that was predicted for q = using bosoniza- 
tion analysis [l8| when changing the scattering lengths 
ao or a,2 (in case of attractively interacting / = 3/2 
fcrmions similar phase transition happens between sin- 
glet and quarteting phases [22Tl2^ ). Our work suggests 
that this phase transition can be induced by external 
magnetic field via quadratic Zeeman coupling, without 
the need to change scattering lengths, thereby present- 
ing certain advantages for experiments. 

Relevant model {or 40 K atoms- To address the nature of 
ground states for the 4-component mixture of 40 K atoms 



we have to break the rotational symmetry associated 
with the hyperfine spin / = 3/2. We relabel the states 
from the / = 9/2 manifold that were used in experiment 
[H as follows, to/ = 9/2 -> 1/2, m f = 1/2 -> -1/2, 
771/ = 7/2 -> 3/2, and 777,/ = 3/2 -> -3/2, so that the 
lattice Hamiltonian is again given by Eq.([T|), but with 
the modified quartic part, 



H 



int,j 



a<0 



(3) 



All interaction constants, U a ^, are different combina- 
tions of scattering lengths af , 0,4 , a% , and aff that we 
take from the table [l[ and spin-changing interaction am- 



plitude U 1 ~ ^_(af 



For q = and t — > corre- 



sponding to infinitely deep lattice, starting from the t = 
state and performing degenerate perturbation theory, one 
can show that a classical (Ising) Neel state is energetically 
favorable. The Neel state breaks translational symmetry 
with the original m/ = 1/2 and m/ = 3/2 components 
paired on odd sites and m/ = 9/2 and to/ = 7/2 paired 
on even sites (or vice versa) . The pairing of to/ = 1/2 
and to/ = 3/2 components minimizes the interaction en- 
ergy per lattice site, and due to constraint: N±/ 2 = Nt 



v 9/2 



and N 3 / 2 = N 7 / 2 other L/2 sites must be occupied by 
to/ = 9/2 and to/ = 7/2 pairs. Activating t introduces 
Neel order due to derealization energy gain. 

To study effects of quadratic Zeeman coupling we use 
infinite system size DMRG (2f| calculations, and moni- 
tor different order parameters depicted in Fig.Q. Ex- 
plicit expressions of non-local parity and string orders 
can be borrowed from the corresponding orders of S=l 
spin chain [26j by identifying Neel up ( mt = 1/2 and 
777/ = 3/2 pairs) and down (to/ = 9/2 and to/ = 7/2 
pairs) components with S z = +1 and S z = —1 states, 
and on-site singlets (initial pairs of BI composed from 
777/ = 9/2 and to/ = 1/2 atoms) with S z — 0. The se- 
quences of phase transitions in strong coupling induced 
by q are shown in FiglU For deep lattices the emer- 
gence of Neel and Haldane phases in the presence of 
quadratic Zeeman coupling is a generic feature of 4- 
componcnt fcrmions with spin-changing processess al- 
lowed for / 7^ 3/2. Haldane phase can intuitively be 
understood as a Neel order diluted with the defects that 
represent on-site singlets of the BI state. The defects 
are created on even number of adjacent sites due to the 
quantum fluctuations on top of the underlying Neel state. 

For shallow lattices, as depicted in FigJTJ a dimer state 
wins. This is consistent with weak-coupling bosonization 
analysis. For the scattering lengths taken from the exper- 
iment in [l[ , the transition from dimer to Neel state hap- 
pens in the strong coupling regime with increasing lattice 
depth, where bosonization cannot be trusted. However 
one can induce similar phase transition in weak-coupling 
by changing interaction coefficients in Eq.©. This allows 
to understand dimer to Neel transition qualitatively as a 
Gaussian transition in the <pt 2 sector. In dimer phase the 
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FIG. 4: (Color online) Quadratic Zeeman coupling (measured 
in units of t) dependence of the order parameters along a cut 
of Fig[T]for U/t = 20. 

fields are pinned as follows: (<p v ) — {4>ti) = {4>t2) = 0, 
whereas in Neel phase: (<p v ) = (</>ti) = and (</> t2 ) = 
[l8j]. With increasing q an Ising phase transition 
happens from Neel to Haldane phase in the 4> v sector 
also captured by the effective model of Eq.((2). The ex- 
pectation values of bosonic fields in Haldane phase are: 
(6 V ) = (4>ti) = and (4>t2) = an d with further in- 

creasing q the singlet/BI state should be recovered, with: 
(0 V ) = ((f>ti) = {4>t2) = 0. Hence, occurence of Haldane 
phase between Neel and BI is natural from bosonization. 
Between Haldane phase and BI (like between dimer and 
Neel) a Gaussian criticality from gapless 0t 2 sector is ex- 
pected, characterized by non-universal exponents. This 
is in agreement with DMRG results failing to capture uni- 
versal exponents at these phase transitions and confirm- 



ing their smooth nature (second order). The complete 
numerical ground state diagram is presented in FigJTJ 

For A^-component fermions at half-filling, with at- 
tractive interactions and 5(7(2) spin rotation symmetry, 
the emergence of charge Haldane insulating phases was 
shown in Q and spin Haldane phase was discussed for 
alkaline-earth atoms in (To| . 

In summary- a rich ground state phase diagram of half- 
filled 4-component repulsively interacting alkali fermions 
can be explored by changing the lattice depth and the 
quadratic Zccman coupling, without the need to change 
the scattering lengths. For the case of / = 3/2 fermions 
dimer and singlet phases are realized, whereas when 4 
components are from / > 3/2 multiplet (like 40 K atoms), 
four distinct gapped MI phases are expected. In partic- 
ular a topological Haldane phase can be stabilized start- 
ing from the BI initial state (with vanishing entropy per 
particle) and allowing spin changing collisions by adi- 
abatically changing the magnetic field. We hope that 
our work will motivate experiments on spinor quantum 
lattice gases to observe spin order beyond the super- 
exchange scale. These MI phases may be created and 
detected within the current state-of-the-art techniques in 
optical lattices. Due to a shallow harmonic trap along the 
lattice the MI phases will occupy the central region as in 
the case of BI [l| . The different phases may be detected 
by different means including monitoring the chirality in 
time of flight [l[ and using Faraday rotation as proposed 
in [I?} • The string-order of the Haldane phase may be 
studied using similar site-resolved measurements such as 
those recently reported for the measurement of nonlocal 
parity order in MI [28| . 
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Quantum Engineering and Space- Time Research) and 
DFG Research Training Group (Graduiertenkolleg) 1729. 
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